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^^ ' We provide a novel mechanism that resolves the Big Bang Singularity present in FRW space- 

^Z^ ' times without the need for ghost fields. Building on the fact that a four-fermion interaction arises in 

Cn , General Relativity when fermions are covariantly coupled, we show that at early times the decrease 

in scale factor enhances the correlation between pairs of fermions. This enhancement leads to a 
BCS-like condensation of the fermions and opens a gap dynamically driving the Hubble parameter 
H to zero and results in a non-singular bounce, at least in some special cases. 

fT^ ' PACS numbers: 

-I— > 
Qj- I. INTRODUCTION AND MOTIVATION 



As is well known the Standard Big Bang cosmology (SBB) successfully predicts the observed large scale expansion, 
thermal properties, nucleosynthesis and the cosmic microwave background in our universe. However, the theorems of 
Hawking and Penrose prove that a curvature singularity exists in the SBB at the "birth" of the universe's space-time, 
yj^ ' Cosmologists have had long held expectations that at the Planck time quantum gravitational effects might resolve 
^^ I the singularity and provide a quantum bridge that connects a collapsing phase to an expanding one; otherwise known 
^^ ' as a bouncing cosmology. The implementation of the bouncing scenario in quantum theories of gravity are still at 
the toy model stage, since a complete description of quantum gravity is lacking. Given this fact, we will argue that 
f~^ ' a non-singular cosmology is possible from the effects of fermions on space-time when the scale factor starts becoming 
00 ', small. 

Most modern approaches towards a unified quantum theory of gravity, such as supergravity [l|, or LQG (3 has found 
the first order formalism of gravity, where one treats the vielbein and connections as independent variables, to be 
the natural starting point, ft is also well known that while passing from the first to the second order formalism (by 
integrating out the torsion field) the covariant coupling of fermions to gravity yields a four fermion interaction in these 
d theories [H, Slj, y, ly, UJ ^ . Now, the BCS theory of superconductivity states that if there exists a weak attractive four 
fermion interaction and a fermi surface, fermions will condense into a new ground state with lower energy comprising 
of spin-0 Cooper Pairs. In the context of a homogeneous and isotropic cosmological space-time^, we will show that 
fermions indeed form Cooper pairs. Moreover, the FLRW scale factor (or the "volume" of the universe) plays the 
role as an effective coupling constant such that at early times the correlations between the fermions get stronger. As 
a result, we show that as we approach the singularity an energy gap opens up, contributing negative energy which, 
as we shall explain shortly, is one of the fundamental requirements for obtaining realistic bounces in a spatially flat 
or open universe. We were also able to explicitly show that at least when the usual (positive) matter energy density 
has an equation of state parameter, to < 1/3, the negative gap energy can drive the expansion rate of the universe 
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^ See however Q for a counter example where even in the presence of the Immirzi parameter, the fermions don't change the classical 

equations. 
^ The interplay between inhomogeneities in the Fermi-gas and the metric, and especially it's implications for CMB remains an important 

open question. 



to zero, leading to a bounce between the contracting and expanding phase. For this mechanism to work, an a priori 
assumption that is required is the existence of a finite non-zero density of the Dirac fermions {i.e. we assume a 
particle/anti-particle asymmetry from the very beginning). We employ robust methods of effective field theory to 
demonstrate the non-singular behavior of this cosmology. We note in passing that similar ideas have been previously 
employed to resolve the black hole singularity and information loss problems 9] . 

In contrast, previous attempts to resolve Big Bang singularity (BBS) has mostly relied on introducing pathological 
ghost fields flO, HH II4I which either violate quantum unitarity and/or leads to catastrophic instability [ij] . In some 
other bouncing and cyclic universe scenarios [ij, [l^, [la], one leaves the singularity unresolved in the hope that 
"quantum gravity" effects will eventually smooth out the singular Big Crunch/Bang transition. Although rigorous 
but involved no-go theorems exist in the literature jl7| as to why it is so challenging to avoid the singularity in General 
Relativity (GR), in the context of a flat'^ homogeneous isotropic cosmology it is easy to qualitatively see this. One 
just has to look at the Hubble equation 

which governs the evolution of the scale-factor a{t) appearing in the Friedman-Lemaitre- Robertson- Walker (FLRW) 
metric 

ds^ ^~dt^ + a'^{t)dxl . (1.2) 

Here, i labels the different components of matter present in the universe. Now, in the bouncing universe construction 
the scale factor goes through a minimum, where d, or equivalently the Hubble expansion rate, H , must vanish. This 
means that some of the matter components must have negative energy, a property hard to realize in conventional 
matter. Ghost fields have negative kinetic energies and can mediate a bounce and therefore have been used in previous 
literature. One may wonder whether ordinary scalar fields with negative potential energy can work? Unfortunately, 
having a negative energy component is not sufficient for a bounce, the negative energy component has to precisely 
cancel the positive matter component at the bounce and then redshift away faster than it, to leave the total energy 
density positive. In other words one requires violation of not only the weak energy condition (p < 0) but also the null 
energy condition {p-\-p <Q) [l7|. As a consequence the energy density has to increase with expansion of the universe! 
It is easy to check that ordinary scalar fields cannot achieve this. 

The reason why an attractive four fermion coupling can give rise to a consistent bouncing scenario is because, firstly 
the interaction energy between fermions gives rise to a negative contribution to the energy density (binding energy 
between Cooper pairs). Secondly, this energy density depends on the "Gap" which in turn, depends rather nontrivially 
on the scale factor (or volume of the universe) via the chemical potential. It turns out that this rather non-trivial 
volume dependence can violate the null energy condition temporarily as needed to have a bounce. On hind sight the 
fact that fermions can violate the energy conditions probably should not come as a real surprise, even "classically" 
fermionic condensates are known to violate energy conditions and give rise to bounces [20|. 

We should point out that in the context of effective field theory there have been some successful efforts in constructing 
toy models involving higher derivative actions which can evade the pro blem of ghosts (for instance involving non- 
local modifications of gravity [2l|, or ghost condensation mechanism [22|) but preserve some of the useful ghost-like 
properties of higher derivative theories to resolve the BBS. (For a more detailed review of various bouncing scenarios 
the readers are referred to [23|.) However, the mechanism that we propose is not a toy model, and has a much more 
generic scope as it essentially only relies on having covariant coupling of fermions to gravity and the presence of a 
Fermi surface. 

A third requirement to achieve a bounce intriguingly, turns out to be the presence of the Immirzi parameter (associated 
with Hoist's generalization of the Einstein-Hilbert action) which also seems to play a key role in the nonperturbative 
quantization of gravity 2]. For instance, in Loop Quantum Gravity (LQG), the quantization of the area operator, 
and the existence of bouncing cosmologies [2J] depend crucially on the Immirzi parameter. What is also suggestive 



^ Although our analysis is not specific to a flat universe and readily applies to open or closed universes, for simplicity and phenomenological 
reasons (WMAP and other measurements strongly constrain the flatness [l8l |) we specialize to the flat case. However, we point out that 
in a closed universe setting the presence of the spatial curvature itself can allow for a resolution of the big bang singularity [191 , but 
obviously such a resolution mechanism will not work for a general manifold, such as an open or a flat universe which makes it somewhat 
less appealing. 



is that like in LQG, we have to rely on a completely nonperturbative effect, the formation of a BCS gap to evade the 
Big Bang. 

This paper is organized as follows: In section II we present the first order formalism of GR and derive the four 
fermion interaction. In section III we develop the cosmological BCS theory and derive the effective potential for 
the cosmological energy gap. In section IV we analyze a special case in details to show how and when singularity 
resolution can occur. Finally we conclude and provide future outlook in section V. 



II. FOUR FERMION INTERACTION FROM GR 

The goal of this section is to summarize the first order formalism of GR in the presence of fermions and derive the 
four- fermion interaction necessary for BCS condensation [5|, |6|. Let us consider a 4-dimensional manifold M. and 
introduce two independent fields: the tetrad, e^, an orthonormal coordinate basis for each point on the manifold, and 
a spin connection A„i/"^ which connects (parallel transports) the tangent spaces at different points of the manifold. 
Note that small-case Latin letters starting from m,n . . . denote spacetime indices, while capital Latin letters starting 
from I, J . . . denote internal Lorentz indices. One can then associate a 4-dimensional metric gmn via 



J „^„ 



(II.3) 



where the Minkowski metric is to be viewed as the metric of the internal space. Internal indices are raised and 
lowered with the Minkowski metric, while spacetime indices are raised and lowered with the spacetime metric. Using 
the connections one can define the covariant derivatives of mixed tensors via usual rules such as 

D„iknl = dniknl + Amn^kpi + v4„i/"'A:„j, whcrC Amn^ = el^e^jA„i/ (II. 4) 

Dm is known as the generalized covariant derivative. As is evident, we use the tetrad and its inverse, e^j to convert 
internal indices to space-time indices and vice-versa. The requirement that the spin connection be torsion free is 
simply A[m„f = 0. 

At this juncture we should briefly discuss the relation between internal and spacetime quantities. Riemannian fields, 
like the metric or the Einstein tensor, live on some substructure of the base manifold M, such as the tangent or 
cotangent spaces, and thus, have a fixed dimension (that of M). Alternatively, gauge fields live in an internal vector 
space, which is independent of the base manifold and could in principle be infinite dimensional. A fiber bundle is 
then simply the union of the base manifold and the internal space, where each fiber is a copy of the vector space 
corresponding to a particular space-time point. Such fiber bundles have an associated structure group or Lie group, 
which qualitatively glues all fibers together. In the context of the first order formalism it is convenient to view m,n . . . 
as associated with the cotangent/tangent spaces, while I, J . . . with internal vector indices associated with the Lie 
Group SO{3,l). 

With the generalized covariant derivative and spin connections we can now define the generalized curvature tensors. 
This is done through the failure of commutativity of the generalized covariant derivatives. One defines 

(II.5a) 
(II.5b) 

where the commutator is short-hand for 

(n.6) 

Note that if the connection is metric compatible and torsion- free (ie, if it is the Christoffel connection), then the 
curvature tensor is simply the Riemann tensor. 

Let us now rewrite the Einstein-Hilbert action in terms of these new variables. Note, however, that we wish to 
work with the trace of the generalized curvature tensor, and not the Ricci scalar, since these two quantities are not 
necessarily equivalent. The Einstein-Hilbert action is given by the well-known expression 

j\,r2 p i\j2 p 

Se = ^ j d^x V^R ^^Jd^x eeTe'jP^-'KLR^n'"'', with P^'^kl = S^^S^^ (II.7) 
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= dlrrAn]/ + [A,n,An]j-^, 
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= d[mAn]p'' + [A,n,Ar^]p'', 


\AjYii A^\j ' 


•= Ami AnK ~ AnI AmK 



where E simply stands for "Einstcin-Hilbcrt" . The second relation can be derived using the identity 

r^Ti en 



R = bi;s",R^j'^, (II., 



We note in passing that in the absence of ferniions, the field equations for the connection gives rise to the compatibihty 
(or zero torsion) condition: 

Di^e^] = (II.9) 

which means that Ami"^ are nothing but the Christoffel symbols determined in terms of the metric and one recovers 
Einstein's GR. This statement remains true even when one includes Hoist's modification to (JII.7P : 

Pk'KL ^ pI/kl ^ S^kS'i] - ^e'-^KL (11.10) 

One can check that once one imposes the compatibility condition (|II.9p . the additional term in Hoist's action, Sh 
(which has Ph rather that Pe in (|II.7p ). vanishes and therefore does not effect the classical equations of motion. 

However, things dramatically change when one includes the covariant coupling of free fermions to GR in the presence 
of the Hoist term: 

S = Sh + Sd (Il.lla) 

Sd = -^ fd^xV^{iWeTT^„.t^ + c.c.), (Il.llb) 

Sd corresponds to the Dirac action with massless fermions '^, where c.c. stands for complex conjugation, -0 is a Dirac 
spinor, and 7^ are 4x4 gamma matrices defined via 

-n= ( ^ ^q] with ai = (1, -a) and ct/ = (1, ct) (11.12) 

where a are just the usual 2 by 2 sigma matrices. With these definitions it is easy to check that the gamma matrices 
obey the following anticommuting algebra 

{7/,7j} = -2r;,j (11.13) 

In pi.llb|) we have also defined 

V; = ^^70 = V'^ f 5 I) (11-14) 

Note that a tetrad based formalism is essential for the inclusion of fermions in the theory, since Dirac spinors live 
naturally in SU{2). Therefore, covariant derivative associated with the Dirac action are not the usual 50(3,1) 
covariant derivatives, but instead are given by 'Dm'4' ■= dm'4' ~ (l/4)^m 7/7j V'- 

Let us now find the structure equations of the fermion-extended theory. In general one can break up the connection 
into symmetric and antisymmetric pieces: 

Ai^^uj'J + C^' (11.15) 

where uj^^ is the torsion-free spin connection satisfying the compatibility condition pi.9p and Cj^^ is the so called 
"contorsion" tensor. The idea is to integrate out the the contorsion tensor which then will lead us to the more familiar 
second order formulation of gravity where the connection are just the metric dependent Christoffel symbols. This can 
be achieved simply by imposing the structure equations obtained by varying the action with respect to the connection. 
Following 1^] and using the identity 

^i^U^K] ^ _^^IJKL^^^^ ^ 2ry^[^7^'l, (H.16) 



The inclusion of a mass term does not affect the conclusions of this paper, but we leave it out to avoid cluttering. 



we can express the contorsion tensor in terms of the axial fermion current, Jg :— ^757^1/) 

eTCmjK = 47rG'-j— - ( -euklJ^ ViljJsK] j (11-17) 

From the above expression for the contorsion tensor it is clear that CmjK is a non-propagating field, its field equations 
do not have any derivatives on it. Thus "integrating it out" is not only equivalent to reinserting its expression (|II.f7[) 
in the full action classically, but also quantum mechanically. Thus the four-fermion contact interaction term that we 
arc going to generate in going from the first to second order formalism is quantum mechanically an exact result. This 
is a key difference from the "effective" contact interaction that one obtains in non-abelian gauge theory where the 
mediating gauge fields do indeed propagate, and therefore the contact term is only a low energy approximation. 

Substituting pi.l7p in Hoist's action we find that the action can be written as Q 

S=SEH+SD[i0]+S^nt. (11.18) 

The first and the second terms are the standard Einstein-Hilbert and Dirac action involving the Christoffel connections. 
However, crucially one obtains a third interaction term, given by^ 



S^nt = ^-kG ( ^^YT J j "" '^"'^^ = j d'^xe -^ (H.19) 

Such four-fermion interactions were already observed in Einstein-Cartan theory (7^ —> 00 limit in pi.l9[) ). although 
they are suppressed by a power of Newton's constant (a factor of 1/k here). In order to form a condensate of fermions 
the Planck suppression will have to be transcended, and this happens when the universe contracts to Planck densities, 
as we shall discover in the next section. 



III. COSMOLOGICAL BCS THEORY 

In this section the starting point of our discussion is the General Relativistic action pi.lSp derived in the previous 
section, but we want to study it in a condensed matter framework, i.e. not in vacuum but in the presence of a gas® 
of fermions. In this case, one has to add the contribution of the chemical potential (/i) to the effective action, which 
at zero temperature reads [2a, |23] : 



S= d 



-j4„ 



xe 



/Vf2 7 7/ 



(IH.20) 



We note in passing that although here we consider four fermion coupling arising from covariant coupling of fermions 
to gravity via the torsion constraint, an "effective" gauge mediated contact interaction term [26j can also lead to 
similar mechanisms, and could be interesting to pursue in the future. 

Our aim is to study in details the possibility of a cosmological BCS-like condensation of the fermions first suggested 
in [28] and to consider its implications towards the resolution of the Big Bang singularity. As discussed in the 
introduction, to have a resolution of BBS, we have to ensure that the Hubble expansion rate can vanish. As a first 
approximation then it becomes sufficient to look at the condensation mechanism on a fiat Minkowski background 
i.e. ignoring the space-time dynamics; one can readily convince oneself that corrections to due to expansion of the 
universe are expected to be 0{H/ jjl) and therefore are suppressed near a possible bounce when H ^ Q. 

Let us therefore focus on the Minkowski space-time action: 



S= d^ 



'iip-f'^dmip - M'07oV' ■ 



T jm 

M2 



(IH.21) 



The sign difference in front of the four-fermion term as compared to what was derived in [^ [QI is simply because our metric has the 
opposite signature. 

Perhaps it is better to call the system as a Fermi-liquid |25| rather than a Fermi gas, as they are interacting. Also, in principle there 
could be several species of fermions, but for simplicity we will only consider one. 



Technically it is simpler to work in terms of two component Weyl fermions, Cf,a, where the left handed fermion sector 
is denoted by F (it contains the left handed fermion and right handed anti-fermion) and the left handed anti-fermion 
sector is denoted by A (it contains the left handed anti-fermion and right handed fermion). Let us first look at the 
free part of the action: 

s^Jd^x [z(c^t^™a„c'' + C^^^™a™c^) - A^(c^k'' - C^^C^)] (111.22) 

A simple and physically transparent way to understand the condensation mechanism is to introduce auxiliary scalar 
(gap) fields, which are proportional to the fermionic bilinears. The so called "gap equation" is then derived by 
integrating out the fundamental fermionic degrees of freedom. As is well-known J26l . |27|, in these integrals the 
contribution of the gap is inversely proportional to the energy of the fermions. Now each of the ^'s describe two 
physics degrees of freedom, corresponding to the positive and negative helicity states —id. a C — a.p ( = ±|p| (. Thus 
we see that the positive helicity state for (-^ and the negative helicity state for C"* {i.e. the left and right handed 
fermions) correspond to energies E = ±||p| — /^| (the positive and the negative sign corresponds to "particle" or "hole" 
like excitations around the Fermi surface respectively) . Thus close to the Fermi surface the energy of these fermionic 
states can vanish and thereby contribute significantly to the gap energy. On the other hand, the antifermionic states 
corresponding to negative helicity and positive helicity states for C^ and ^ respectively, always have non- vanishing 
energies ^ |p| + /i and hence their contribution towards the gap is suppressed (see for instance [26] for a discussion on 
this point). In the following analysis therefore we are only going to focus our attention on the fermionic states and 
couplings between them. 

A simple way to project out the anti-particle states, as advocated in [26|, is to first go into the momentum basis: 

C(^) = felCpe'P" where Cp ^ J d^x C(x)e~'P^ . (III.23) 

Then for a given momentum mode choose a reference frame where the momentum is aligned along the positive z-axis 
to compute the different terms in the action, and covariantize the final result. It is easy to see that the free action 
becomes 

j d'xc,.. = / (01 [^p ^(^° + ^p)C + C^(P° - ^p)C] ^iie^e ^p ^\p\-f^ (ni.24) 

where C,^ and C,^ are now single component fermions corresponding to the left and right handed fermionic states. 
When pis along the direction of the positive (negative) z-axis they are the upper (lower) and lower (upper) components 
of Q^ and C,^ respectively. 

The four fermion term simplifies as the following: Firstly 

J'" Jf„ = V^nie^-a-^e + C^ta'"C^)(C^t^"C^ + C^^a"C^) = 27,„„C^t^™^F^At^n^A 

The second equality stems from the fact that we are only interested in coupling between C,^ and C, , or more precisely 
between C^ and C,^ (all other terms when projected onto fermionic states will be of the form (C^C^^)^ or (C'^^C^^)^ 
which must vanish due to anti-commutavity). This further enables us to consider only the diagonal a™ matrices (the 
off-diagonal ct™'s only contain couplings between fermionic and anti- fermionic states): 

Having found the desired left-right coupling, one can now introduce the auxiliary fields (A, A*), so that we can rewrite 
the four fermion coupling as 



,.^(C-'0(C-'C')' , i^^ 



2 






(III.25) 



It is clear that the auxiliary field, A ~ C^^C^; ^md a non-zero value for it would signal a cosmological BCS-like 
condensation! In order to find such a nontrivial value for A, one can now takes recourse to a mean-field approximation 
where A is treated as a constant "gap" . 

To obtain the "gap equation" (or the effective theory of A) we have to integrate out the fundamental degrees of 
freedom, the fermions. In the momentum basis the fermionic action looks like 

^d'xC,^^ = / (04 {(i'^v' + ep)C^ + C^Hp" ~ e,)C^' + A^^^C^ + AC^t^^] (^1.26) 



Or, in matrix notation 

d^x /:,er = / (04 iCp'^Cp') ( ''"a^'^ p^t,^){lf)^Jd'p {C^KC^')A, ( I ) (III.27) 

As usual in the path integral when one integrates over all the fermionic variables ^p's, one ends up with a fermionic 
determinant. In general for a fermionic path integral we have (see |26l. |27|, for instance) 



Z = / P^e'* = / V^ exp 



A e^MpC 



pSp 



e^pJ^tr{]nMp) (111.28) 



where ^ in general can be a A^ component spinor and Mp is an N by N matrix. The effective action Sett (in Quantum 
Field Theory literature this is often reffered to as Fef f ) is defined in terms of the path integral via 

Z = e'^^" (III. 29) 

Thus for our path- integral we obtain the following effective Lagrangian for A (for a more detailed derivation see [23| ) 

/:eH(A) = -^J ^ trilnAp) — AA* (III.30) 

The above effective action can actually be evaluated exactly. Focussing on the first (non-perturbative) term which 
was obtained by integrating out the fermions we find 

A 



^^P / 2 , A? f '^^P i^P + A*) 



/:non(A) ^J J^^/4 + ^'- J ^ ' v.. "^ [^-l + A^] (111.31) 

where in the last step we have performed the angular integrations and we are considering modes (electrons or holes 
with positive or negative ep respectively) around /i, A is the UV cutoff used to regulate the integrals. The above integral 
is formally divergent. However, employing the above cut-off regularization scheme we first of all can subtract the 
quartic divergence by requiring >Cnon(A) to vanish when A = 0, or in other words just subtract the A = piece from 
£non(A)^ Unfortunately, a logarithmic divergence remains in the form of an undetermined mass parameter e^ (= 2A) 
due to the presence of the four fermion interaction, an artifact of the intrinsically non-renormalizable torsion-gravity 
theory. Even on dimensional grounds it is clear that such a cut-off must be present. In condensed matter systems, 
this cut-off corresponds to the Debye frequency related to the lattice spacing. It seems likely that in case of gravity 
Cr would similarly be related to the fundamental quantum (discretization) of space-time as for instance observed in 
Loop Quantum gravity [3]. 

Putting everything together therefore we have the final effective potential given by 

M3A2 1 



^tot =->Ceff(A) = 



A4 ^2^2 /^2 X ^ 

2.2 i6+V+^ --^ r^^z 



(III.32) 



This has one arbitrary parameter e^ as promised, along with the inverse coupling constant Air- It is clear that the 
above potential has a minimum at say A — Aq given by 

c^V^ot 9 A2 /A2-2u2\ Ao 



^ The quartically divergent piece is nothing but the contribution of the fermions to the vacuum energy. This is a specific illustration 
of the famous "cosmological constant problem" : Quantum loop contributions to the cosmological constant are generally known to go 
as ~ 0{Mp) for a theory with a Planck cut-off scale, but observations suggest that we are living in an universe with a much smaller 
cosmological constant ~ 10~^-^'^Mp. This implies an incredible cancelation of a hundred and twenty orders of magnitude between the 
different contributions to the vacuum energy. Since in this paper we are not trying to address the cosmological constant problem, we 
take the usual approach and assume that some unknown mechanism is indeed responsible for such a cancelation/ "renormalization" of 
the vacuum energy down to it's extremely small observed value, ~ 0{mev'^). Since the bounce in our model is governed by much higher 
energy scales ~ Mp, such a tiny cosmological constant will play no role. If on the other hand the vacuum energy at the bounce was 
large and only became small later, it can modify the nature of the bounce. However as we shall see later, whether we have a bounce 
or not only depends on the equation of state of matter, ui, and in particular as long as oj < 1/3 the Big Crunch/Bang singularity is 
resolved via the bounce. Since the cosmological constant has an equation of state uj = —1, it's existence does not pose a problem for 
our mechanism. 



One can think of the above equation as specifying Ap — Ao(/x). It is useful to check that we recover the usual behavior 
of the Gap in the weak coupling BCS limit 26|,|23] when the fermion gas is dilute. For A ^ /i, e,., (|III.33|) tells us 



A « e^ exp 



-Af^TT^ 



2/i2 



(III.34) 



the familiar exponentially suppression of the gap appears. 
Returning now to the general discussion, we can also calculate 



Kin = V;ot(Ao) 



A2 r 

47r2 



3A2 A2 Ao 
■ H In 



M 



(III.35) 



However, the potential energy that we have calculated includes the contribution from the chemical potential as well. 
This is the "extra" /i dependent term in the Lagrangian (|III.21|) : note that the number density, n, of fermions is given 
by [13 






dV^, 



dfi 



(III. 36) 



In order to do cosmology what we really need is the total energy density which is therefore related to Vtot via 

= ^tot + M" (111.37) 



In particular we need to know the dependence of the gap energy density as a function of the scale factor/number 
density. This is provided by the implicit relation between the number density and the chemical potential pil.36[) : 



ov... ^^l^,_^,^ 



dfi 



Thus putting everything together we have 



= V^tot - Ai- 



ai4c 



dfi 



47r2 



27r2 



3Ag 



M 



Ao 



M_v)...f 



(III.38) 



(III. 39) 



Using pil.33|l . (jIII.38|) and (JIII.39P one in principle knows the dependence of the energy density as a function of the 
fermion number density, pgap = Pgap(«-)- 



IV. BEC AND BOUNCE 



For the present paper, we are interested in homogenous and isotropic FLRW type cosmology focussing on the regime 
H <^ fi, relevant for discussing nonsingular cosmological models. In particular what we would like to demonstrate is 
that in the presence of a fermion condensate, the universe does not contract all the way to a singularity where the 
energy densities go to infinity, but rather bounces back at some finite energy density to an expanding phase thereby 
resolving the singularity. Now, in a realistic universe the matter content of the universe will be varied, but for the 
purpose of illustration here we will assume that the matter content of the universe, apart from the condensate, is 
described by an ideal fluid, pmat, with an equation of state parameter to: 



P.at = Poa-"^'+'^^ 



(IV.40) 



To check the robustness of the bounce mechanism we should check whether or not the universe can bounce for any 
equation of state parameter obeying the weak energy condition, i.e. —I < uj < I. In other words, we want to see 
when the Hubble equation 



H^ = 



1 



3M2 



(Pgap + Pmat ) 



Ptot 

3M| 



(IV.41) 



allows for a bouncing solution. In order to solve the above equation we need to know how the energy densities evolve 
as a function of the scale factor. For matter, since it is treated as an ideal gas this is given by pV.40p . For Pgap, one 



docs not have an explicit expression in terms of the scale factor, but it is implicitly defined in the following way: First 
of all we know that the number density of the fermions scale as inverse volume: 

n ^ noa~3 (IV.42) 

Now, (IIIL33P implicitly determines Aq ~ Ao(/i), so that pil.38p can be thought of as relating the number density 
and the chemical potential. In other words, piI.SSp and (jIII.38|) let's us determine the chemical potential and the gap 
as a function of the number density. Since the energy density is given in terms of the gap and the chemical potential, 
(IIIL39P implicitly determines the energy density in terms of the number density. Finally using pV.42p we, in principle 
can determine the energy density as a function of the scale factor. 

In order to demonstrate that the universe indeed bounces, we have to show two things: (i) the total energy density 
must vanish, ptot = 0, at some scale factor, say a = ab.(ii) Also, d > at a = a;, to ensure that we indeed have a 
minimum (and not a maximum as happens in a turnaround) for the scale factor. Differentiating the Hubble equation 
we find that the latter condition is equivalent to showing that dp/ da is positive at the bounce point, in violation with 
the null energy condition. In principle, one can scan the entire parameter space^ {po^ui,nQ^er,Mr} numerically to 
determine when the singularity is resolved, but this a rather involved and challenging task which we leave for future. 
In this manuscript, we want to consider a particular limit where the expressions simplify considerably. Let us consider 
the strong coupling BEC (Bose Einstein Condensation) limit [27|, M^ -^ 0, and when A ^ e^- In this case (|III.33|) 
simplifies to give us a rather simple relation between the chemical potential and the gap: 

A^ w 2^2 (IV.43) 

The expressions for the energy and the number density also simplifies considerably: 

A3 A 

n ^ naa-^ « -^ In — ^ Aq - a'^ (IV.44) 

O A 4 A 

and pgap « -^In^-^-a-^ (IV.45) 



87r2 



er 



where in inferring the scale-factor dependence of the gap and the energy density we have ignored the logarithms. 
Thus approximately, the gap energy density behaves as radiation with a negative sign. It is worth noting that Casimir 
energy calculations are also known to lead to similar negative radiation like behavior for massless minimally coupled 
fields [2^ , as the case we are considering. 

What the above behavior suggests is that as long as ui < 1/3, i.e. the matter energy density blue shifts (during 
contraction) slower than radiation, we will have a bounce, see figure IIVI Since ptot is always positive, the matter 
energy density will always dominate over the gap energy density, but during contraction since the gap energy is blue 
shifting faster as compared to matter, it will eventually catch up with matter and precisely cancel it at the bounce 
point. After the bounce, in the expanding phase the gap energy will dilute faster than matter ensuring that ptot 
remains positive. In our case, one finds that for 



Ptot ^^ Pbounce 



, 3(l+w) 

a \ ^ ' a 



we have 



and thus as argued before we satisfy the bounce criteria when uj < 1/3. 



dp 
da 



1 - 3w (IV.46) 



A few remarks are now in order. Firstly we realize that in the strong coupling BEC limit that we are considering 
here, whether we will have a bounce or not only depends on uj. In general however, the dependence will be much 
more complicated and only future numerical exploration will be able to provide a comprehensive answer. Let us next 
look at the robustness of the mechanism in resolving the Big crunch/bang singularity. Usually it is believed that in 
a contracting universe at high energy densities the universe will come to be dominated by radiation which goes as 
a~^ and has an equation of state lu = 1/3. Whether we can have a bounce, or not, in this limiting case of w = 1/3 
crucially depends on the log like behavior of the gap energy. By inspection one finds that because of the presence of 
the logarithm in (|IV.44p . Aq blue shifts slightly slower than I /a, and again by inspection from (|IV.45p one deduces 



Actually there is a gauge redundancy in the set of parameters since one can always re-scale the scale-factor. 
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-0.10 



0.10 



FIG. 1: Plot of the scale factor (red curve) and the gap (green curve) as a function of time as the universe bounces. This plot 
corresponds to tj = 0, and e^, Mr <^ Aq. 



that 



nA 



^A, and therefore blueshifts slightly slower than radiation. This means that unfortunately in 



the presence of pure radiation (if other forms of matter are present, the situation may change), one cannot bounce 
back. Although this is a drawback of the proposed mechanism, according to string theory this may not necessarily be 
a serious problem. According to string theory at high enough energy densities, close to the string scale, the thermal 
state of the universe is no longer describable by massless radiation modes. Rather, one enters a stringy "Hagedorn 
phase" where even the massive modes are excited and are in thermal equilibrium [30|. In this phase, the total matter 
energy density (now consisting of both massless and massive modes) actually blue (red) shifts as a~^, and therefore 
the presence of condensates can indeed produce a bounce. 

One may also be worried about the problem of Mixmaster chaotic behavior near the bounce, as the anisotropics are 
known to grow as a^^ as the universe contracts. Since pgap blueshifts slower than anisotropics, our bounce mechanism 
is indeed susceptible to this problem. However, recent work on cyclic cosmologies involving many bounces have been 
shown to, at least, ameliorate the problem p3l| . 

Finally, what about the maximal energy density at the bounce point, what parameters does it depend on? As remarked 
earlier, since one can arbitrarily re-scale the scale factor, only a combination of the parameters no,po is physical. In 
the special case when a; = (as one expects in the Hagedorn phase) for the limiting regime under consideration, one 
can check that 



Pbounce 



256 pl,(t) _ 256 



Pt 



277r2 



'it) 



277r2 



(IV.47) 



To summarize, in our picture the universe starts out with a dilute gas of fermions with an exponentially suppressed 
gap (jIII.34p . This would be completely overwhelmed by finite temperature effects and we will just have a theory of 
ordinary non-interacting fermions. However, as the universe contracts and the number density n increases, so does the 
gap energy. (Increasing fi is equivalent to decreasing M^, or increasing the coupling of the four fermion interaction. In 
this sense the volume or the scale- factor controls the strength of the interaction.) Eventually at extreme high energy 
densities Cooper pairs are formed (superconducting phase), and the negative interaction energy starts to cancel the 
positive kinetic energy contributions more and more. Finally when A > e^, we expect these Cooper pairs to condense 
more and more in the new ground state and the interaction energy can completely cancel the matter energy density 
leading to the bounce. 
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V. CONCLUSION 

In this paper we have presented a novel physical mechanism which self consistently resolves the initial big bang 
singularity. The crucial ingredient relies on the transient violation of the null energy condition during the bounce. 
The interactions responsible for the negative gap energy at the bounce is enhanced as the universe contracts to the 
singularity. Typically it is difficult to obtain a bounce without the introduction of dangerous ghost states, but in this 
mechanism the negative energy arises from the binding energy associated with the bound state of the Cooper pair. 
Moreover, the gap only becomes significant at early times and consistently redshifts away at late times, in other words 
the late time cosmology is protected from the negative energy contribution that generated the bounce. In a special 
region of the parameter space we were explicitly able to obtain bouncing solutions in the presence of matter provided 
it's equation of state parameter was less than one-third. 

While our results are promising, a few developments are in order. First, we have not included perturbations of the 
condensate which is expected to affect the homogeneity and isotropy both at early and late times; how robust is 
the BCS bounce in the presence of inhomogeneities? Second, what would be the observational consequences of this 
mechanism? How would inhomogeneities responsible for large scale structure be seeded in this model? Since we 
have a bounce, we should be able to employ the techniques and physical picture in close semblance to the Ekpyrotic 
scenario. A plausible structure generating scenario might arise from inhomogeneous excitations of the condensate 
something we plan to investigate in the future. A promising route would be to employ the techniques developed 
in |32l | to propagate perturbations (possibly generated by an ekpyrotic scalar field) across the gap-mediated bounce 
as the gap only depends on the overall volume of the universe and is not affected by fluctuations in the metric, a 
pre-requisite for the success of the mechanism advocated in [32| • 
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